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Abstract 

Using the CLEOIII detector at CESR, we have measured the branching ratio of the de- 
cay H° ^ pK-K--K+ relative to ^ H-7r+. We find ^ pR- K--K+) /B{rP^ 
S~7r"'") = 0.35±0.06(stat)±0.03(syst). In the resonant substructure of this mode, we find ev- 
idence for SO decays to plC^iSmfR-, and measure 6(2° ^ plC* {S<^2f K') ■ B(K* (892^ 
K-'iT+)/B{Bl S-7r+) = 0.14±0.03(stat)±0.01(syst) and B{Bl pK' R- tt+) / Bi^E^ 
E;-7r+)= 0.21±0.04(stat)±0.02(syst) for the non-K*(892)° 2° -> pR-R-TT+ decays. This note 
has the revised numbers on the branching ratios with improved secondary vertex finding algo- 
rithm. 
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In the past decade, singly-charmed baryons (consisting of one heavy quark and two hght 
quarks Qqiq2) have been of interest to many phenomenologists working in the realm of Heavy 
Quark Effective Theory P. The heavy charm quark acts as a heavy nucleus and the light 
di-quark moves around it, analogous to the hydrogen atom. The CLEO experiment has 
discovered many new charmed baryons and measured their properties; in particular it has 
measured many of their relative branching fractions. The study of charmed baryon decays 
is complicated because they can proceed by three distinctly different processes; external 
IV-emission, internal VT-decay, and ly-exchange. Disentangling the contributions of each of 
these processes requires the measurement of as many different decay modes as possible. 

Improvements in the particle identification in the CLEO III detector with the intro- 
duction of the RICH (Ring Imaging CHerenkov) sub-detector [3], have made it possible to 
search for decay modes previously contaminated by huge combinatorial background. This 
paper concentrates on one such decay mode of the S° (the csd charmed baryon, discov- 
ered by CLEO 0]), namely "EP^^pK^ -n^ , and its substructure. To measure the relative 
branching fractions we use s-S~7r^ as the normalizing mode. Charged conjugation is 
implied throughout the text. The only previous observation of the pK~ K~7r~^ final state 
was made in 1990 by the ACCMOR Collaboration who observed four E^^^pK'W (892)° 
decays; but there was no information on their rate. 

The data for this analysis were collected using CLEO HI detector based at the Cornell 
Electron Storage Ring (CESR) taken at and near the T{2S), T(3S'), and T(45') resonances. 
The integrated luminosity corresponds to 7.2 fb~^. In the CLEO HI detector configuration, 
the innermost tracking device is a four-layer double-sided silicon vertex detector surrounding 
the beam pipe. Beyond this vertex detector is the main cylindrical drift chamber 0], with the 
inner 16 layers being axial and the outer 31 layers having a small stereo angle. The tracking 
system is immersed in a 1.5 T solenoidal magnetic field and measures the momentum and 
specific ionization {dE/dx) of charged particles. Outside the drift chamber is the RICH 
sub-detector |2.] consisting of two concentric cylinders. The inner cylinder comprises LiF- 
crystal radiators on a carbon-fiber shell, while the outer one consists of thin multi-wire 
proportional chambers filled with a gas mixture of methane and TEA. Ultra-violet Cherenkov 
photons interact with the TEA emitting electrons that are multiplied by the chamber and 
whose image charge is sensed on pad detecting planes, thus localizing the photon positions. 
Between these two layers is a 16-cm thick gap which allows the Cherenkov cone to expand. 
Surrounding the RICH is the 7800 Csl crystal calorimeter for the identification of photons 
and electrons. Beyond the crystal calorimeter is the superconducting solenoid magnet and 
a muon detector system. The crystal calorimeter and the muon detectors are not used in 
this analysis. 

To identify hadrons, we combine information on the specific ionization (dE/dx) measured 
in the drift chamber and likelihoods obtained from the RICH detector. The RICH likelihood 
is formed for each hypothesis [i = it, K, p) by using the measured positions of photons that 
are located within 3 standard deviations of their expected positions for the measured track 
momentum. The likelihood function for each particle hypothesis, i, is defined as, 

U= X{[G{ef^\eTn + (1) 

where G is a Gaussian-like probability function of observing the photon at an angle 
6'°^* with respect to the expected Cherenkov angle 6^^^^ for particle type i, and S is a flat 
background probability function. The details of RICH identification are given elsewhere ■ 
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Similarly, using the dE/dx information, we construct a quantity for the different hypotheses 
{i = 71, K, p) as Si, which is the difference between the measured and expected dE/dx 
for that hypothesis, expressed in units of its standard deviation. The RICH and dE/dx 
information for each pair of hypotheses is then combined to form functions, for example, 

^X\P - vr) = -21ogLp + 21ogL^ + _ (2) 

In this example real protons peak at negative values of Ax^ip — vr) whereas real pions 
tend to have positive values. If there is no RICH information available, or if the particle's 
momentum is less than 1.0 GeV/c (for a proton), or 0.5 GeV/c (for a kaon), we use only 
dE/dx information to form the Ax^ function. 

All the primary charged tracks are required to have a distance of closest approach to the 
beam position in the r — (p plane of less than 5 mm and of less than 5 cm along the z axis. 
Typical beam dimension is 300 microns, 100 microns, and 10mm in x,y, and z, respectively. 
In the CLEO environment charmed baryons do not have well separated decay vertices. We 
require that the scaled momentum of the charmed baryon candidate, Xp, be greater than 0.5. 

Here Xp = P/ \J El — M^, P and M are the momentum and mass of the candidate, and Ei, is 
the beam energy. This requirement dramatically suppresses combinatorial background, and 
dictates that the observed charmed baryons are the result of continuum production rather 
than from B meson decays. 

Candidates for the S°— >-pii''K^7r+ decay are reconstructed by combining one proton 
candidate, two kaon candidates, and a pion candidate. For the proton identification, we 
require Ax^{p — K) < —4 and Ax'^{p — tt) < —4. Similarly for each kaon identification, we 
require Ax'^{K — p) < —4 and Ax'^{K — tt) < —4. Pions are selected with a loose dE/dx 
criteria jS"^! < 5. All charged tracks are required to have momenta in excess of 100 MeV/c. 
Once the four charged tracks are selected we kinematically constrain them to come from 
a common vertex. The invariant mass distribution of the pK~ K~tt^ candidates is shown 
in Fig. n The data are fit to a Gaussian signal function and a second-order polynomial 
background shape. The fit yields a signal of 148±18 events at a mass consistent with previous 
measurements of the j3] and a fitted width of 4.1±0.5 MeV/c^, consistent with ejmected 
resolution of 4.5 MeV/c^ obtained using a GEANT-based Monte Carlo simulation [9|. 

In this multi-body final state, pK~ K~tt^ , we also search for resonant sub-structure 
-ft"* (892)" by calculating the invariant mass of each of the kaon candidates combined with the 
pion. Figure 121 shows the sideband-subtracted f^'Tr"*" invariant mass using those pK~ K~tt^ 
combinations within 3 standard deviations of the S° mass peak. Candidates for the signal 
(sidebands) are selected within the mass region of 2458.3 - 2483.1 MeV/c^ (2417.6 - 2442.4 
or 2498.3 — 2523.1 MeV/c^), as shown in Fig. ^ The low mass tail in this distribution is 
due to the fact that we combine both kaons with a pion. The correct combination appears 
in the peak region, and the incorrect once forms the low mass tail. This distribution is fit 
to the sum of two shapes. The first one is the ii'*(892)'' signal shape, which is generated 
using a 3-body phase space simulation of s>j9iir~iir*(892)°, with K* {^92)^ ti^ . For 
the second , we use a four-body non-resonant simulation of 'E^^-^pK~ K~ t:^ , since the data 
show no evidence of any other narrow resonances. In particular, we also searched for the 
two-body decay of the H°^A(1520)Z^(892)°, where A(1520)->pfs:- and W{m2f^K-'K+, 
and found no evidence for this mode. Moreover, broad resonances do not yield a statistically 
significant difference in shape than that from non-resonant production. We therefore model 
all the non-i^r*(892)° contributions using this non-resonant production model. We fit this 
plot to the sum of two shapes from Monte Carlo simulation, one obtained for non-i^*(892)'' 
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"E^^^pK' K~ decays (dashed histogram) , and the other for resonant S^— >p/r~i^*(892)° 
decays (dotted histogram), where i^'*(892)°— >i^'~7r"'". In the fit the two normahzations are 
constrained to add to unity. The statistics are too poor to extract any possible interference 
effects in the invariant mass distribution, and is therefore not considered here. The measured 
resonant and non-A'*(892)'' fractions are found to be 0.39±0.06(stat) and 0.61±0.06(stat), 
respectively. Thus, of the total 148±18 fitted S° {—>-pK~ K~7i~^) candidates 58±ll(stat) 
events are contributed by resonant pK^ K*(892)^ decays and 90±14(stat) events are con- 
tributed by non--ft'*(892)° decays. The x^/dof of the fit (shown in Fig.|21) to the resonant and 
non-i^*(892)° components is found to be 44/45, indicating that our fitting function, which 
includes a contribution from only one resonance, is a satisfactory one within the available 
statistics, and we believe that the systematic uncertainty may be neglected. 

The reconstruction efficiency of these two final states are estimated using the Monte Carlo 
simulation. Within the kinematic region Xp > 0.5, we find that both final states have a re- 
construction efficiency of (23±1) %. The corresponding efficiency-corrected measured cross- 
sections times branching fractions (cr ■ B) are 91±12(stat)±8(syst), 37±7(stat)±3(syst), and 
54±9(stat)±5(syst) fb for all pK-K-7r+, F^(892)° resonant (F^(892)°^K-7r+ only) and 
non-ir*(892)*' pK~ K^tt'^ , respectively, where all the measurements refer to that part of the 
momentum spectrum with Xp > 0.5. 

As the production cross-section of S° baryons is unknown, we do not have a measure 
of the absolute branching fraction of any S° mode. Instead, we measure the branching 
ratios of these new modes with respect to that of the well-estabhshed decay S° — > H~7r+, 
where S~— i>A7r~. Both the H~ and its daughter A have long flight paths, with cr values 
of 4.91 and 7.89 cm, respectively. Therefore, in the H°^S^7r+ decay chain we have two 
vertices significantly detached from the beamspot. The A sample is selected by vertexing two 
oppositely charged tracks. The protons from the A decays, which have the higher momentum 
of the two daughters, are required to be consistent with a proton hypothesis {Ax^{p — K) 
and Ax^(p — vr) < 0). Background is further rejected by requiring the daughter tracks from 
the A to be inconsistent with coming from the beam interaction point. The A candidates 
within 5 MeV/c^ (3 standard deviations) of the nominal PDG 8] mass (1115.68 MeV/c^) 
are then kinematically constrained to this mass and combined with an appropriately charged 
track to form the S~ candidate. The A decay vertex is required to be at a greater distance 
from the beamspot than the S~ decay vertex. The S~ candidates are also required to have 
a flight distance of 3 mm or more before decaying. Pions from the S~ baryons are required 
not to come from the interaction point, by requiring the of the fit (if forced to come 
from the interaction point) to be greater than 3. Those S~ candidates within 7 MeV/c^ (3 
standard deviations) of the PDG [8| mass (1321.31 MeV/c^) are kinematically constrained 
to this mass and are used for further analysis. Finally a charged track consistent with the 
pion hypothesis is combined with the H~ candidate to reconstruct the S° candidate. A fit to 
the S'TT"*" invariant mass distribution returns a signal yield of 182 ± 18 and a mass consistent 
with previous measurements as shown in Fig. El The S^^S~7r"'' reconstruction efficiency 
is measured to be (9.8±0.3) %. The measured branching fraction times the cross section 
(a ■ B) for the S-7r+ mode is 260±26(stat)±23(syst) fb for Xp > 0.5. 

The measured relative branching fractions of the S^—^pK^K^ir^ modes are tabulated in 
Table □ 

We investigated several sources of uncertainty, including background shapes, signal width, 
Monte Carlo statistics, charged particle identification, and A and S~ reconstruction. The 
dominant uncertainties arise from the S~7r"'" mode, due to the reconstruction of the displaced 
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TABLE I: Measured branching fractions of the 'E^-^pK'^ ir^ mode relative to that for 
H[?^H^7r+. The errors after the values give the statistical and systematic uncertainties, respec- 
tively. 

Ratio of modes Relative branching fraction 

„ J ^ 0.35 ± 0.06 ±0.03 

B0i^pK~K*{mym*mf ^ K~7r+) Q^^^ ^ 0.03±0.01 



B0^^->pK-K--K+) No K*(892)° 



0.21 ±0.04 ±0.02 



vertices. 

To estimate the uncertainty due to our assumptions in the shape of the background, we 
tried both first and second order polynomials to describe its shape. The fitted yield changed 
by 3%, which we take as our systematic uncertainty from this source. The systematic 
uncertainty due to the imperfect understanding of the signal resolution is taken as the 
difference in signal yield using a floating width and a width fixed to the value found from 
simulation (3% for the branching ratios, 6% and 8% for the pK~ K~ tt^ and the "E'tt^ modes, 
respectively). We assign 4% and 5% systematic uncertainties to the 'EP^—^pK~K^Ti^ and 
H°^H~7r+ decay modes, respectively, due to the finite statistics of the Monte Carlo samples. 
As the number of tracks for the observed (numerator) and the normalizing (denominator) 
modes are the same, the 1% per track uncertainty in basic track-finding cancels in the 
calculation of the branching ratio, as does the uncertainty in the luminosity (2%). Systematic 
uncertainties in the charged particle identification are investigated using samples of protons 
and kaons from the A^pvr" and Ki.-^pK~TT^ modes, respectively. The study demonstrated 
that the uncertainties in the particle identification for protons and kaons in the relevant 
momentum range are 3.2% and 2.4%, respectively. As there are two kaons in the final 
state, we assign 4.8% uncertainty to the kaon identification. Based on a study of displaced 
vertex finding in data and Monte Carlo we assign systematic uncertainties of 6% to the A 
finding and 2.8% to the S~ finding; these numbers include the extra uncertainty in track- 
finding for low momentum, large impact parameter tracks. All systematic uncertainties are 
summarized in Table ITTl The total systematic uncertainty, obtained by adding the individual 
contributions in quadrature, is 9%, and considerably less than the statistical uncertainty. 

In conclusion, we have measured the branching fraction for "E^^^pK' K~ Tt'^ relative to 
>-S~7r"'". We find that 39±6(stat)% of the signal proceeds via the resonance sub-structure 
pK-li*{m2f, where Z*(892)°-^ir^7r+ with the remainder being the non-K*(892)° de- 
cay pK~ K'TT^ . Results are given as a ratio normalized to the H^tt^ rate. The mea- 

A u u- ( B{pK-K-n+) B{pK-W{m2fyB(K^{Smf^K-TT+) , 
sured branchmg ratios for — ^v> ^ — r — ■, — ^ v-, ^ — r ■, and 

B(pK K ^^j^mK*{892) 0.35±0.06(stat)±0.03(syst), 0.14±0.03(stat)±0.01(syst), 

and 0.21±0.04(stat)±0.02(syst), respectively. This is the first measurement of a S° de- 
cay mode where both of the s quarks in the final state are part of mesons. It is possible 
that such four-body decays proceed via external ly-decay, internal W-decaj or Vl^-exchange 
decay diagrams. Resonant decays such as pK* {892)^ K~ , which have no vr"^ in the final 
state, cannot decay via external W-decaj. Their observation is not surprising, as many such 
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TABLE II: Systematic Uncertainties: The systematic errors listed under Ratio are for the relative 
branching fractions and the errors tabulated in the third and fourth column are pK~ K'n'^ and 
S^TT"*" modes, respectively. 



Source 


Uncertainty (%) 




Katio 


pj\ A TT ' 




: : 

Luminosity 




o 
z 


o 
z 


iracK Jxeconst ruction 




A 


A 


R £1 cAc &VCW 1 n H Q n q t» *i 


3 


3 


3 


Signal Width (Monte Carlo) 


3 


6 


8 


pK~ K~TT'^ (Monte Carlo statistics) 


4 


4 




'E~n~^ (Monte Carlo statistics) 


5 




5 


Proton ID 


3.2 


3.2 


3.2 


Kaon ID (two kaons) 


4.8 


4.8 




A reconstruction 


6 




6 


H~ reconstruction 


2.8 




2.8 


Total 


9.8 


9 


9 



A+ and decays modes have been discovered, and is a further indication that external 
VT-dccay diagrams do not dominate in charmed baryon decays. 
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FIG. 1: Invariant mass distribution of pK K vr^ candidates in CLEOIII data. The fit to the 
above mass distribution yields 148±18 signal events. The signal band (2458.3 - 2483.1 MeV/c^) 
is defined within the solid hues and the low (2417.6 - 2442.4 MeV/c^) and high ( 2498.3 - 2523.1 
MeV/c^) side bands are defined by the dashed lines. 
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FIG. 2: The K^tt^ invariant mass in CLEOIII data for H[! candidates within 3 standard deviations 
of the PDG value. The two sideband contributions have been subtracted. Dots with error bars are 
data points, the dashed histogram is the non-K*(892)° contribution and the dotted histogram is 
the resonant contribution. The sohd hne histogram is the sum of the two contributions. 
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FIG. 3: Invariant mass distribution of H vr^ candidates in CLEOIII data. The fit to the mass 
distribution yields 182ibl8 signal events and the fitted mean is consistent with the nominal 
PDG mass 0. 
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